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Abstract 



Electromagnetic modeling provides an interesting context to present a link 
between physical phenomena and homology and cohomology theories. Over the 
past twenty-five years, a considerable effort has been invested by the computa- 
tional electromagnetics community to develop fast and general techniques for po- 
tential design. When magneto-quasi-static discrete formulations based on mag- 
netic scalar potential are employed in problems which involve conductive regions 
with holes, cuts are needed to make the boundary value problem well defined. 
While an intimate connection with homology theory has been quickly recognized, 
heuristic definitions of cuts are surprisingly still dominant in the literature. 

The aim of this paper is first to survey several definitions of cuts together with 
their shortcomings. Then, cuts are defined as generators of the first cohomology 
group over integers of a finite CW-complex. This provably general definition has 
also the virtue of providing an automatic, general and efficient algorithm for the 
computation of cuts. Some counter-examples show that heuristic definitions of 
cuts should be abandoned. The use of cohomology theory is not an option but 
the invaluable tool expressly needed to solve this problem. 
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1 Introduction 



There is a remarkable interest in the efficient numerical solution of large-scale three- 
dimensional electromagnetic problems by Computer-Aided Engineering (CAE) soft- 
wares which enables a rapid and cheap design of practical devices together with their 
optimization. 

Electromagnetic phenomena are governed by Maxwell's laws [T] and constitutive 
relations of materials. This paper focuses on the numerical solution of magneto-quasi- 
static Boundary Value Problems (BVP) — also called eddy-current problems — which ne- 
glect the displacement current in the Ampere-Maxwell's equation [1], [2], [3]. This well- 
studied class has quite a big number of industrial applications such as non-destructive 
testing, electromagnetic breaking, metal separation in waste, induction heating, metal 
detectors, medical imaging and hyperthermia cancer treatment. 

The range of CAE applications is sometimes bounded by the high computational 
cost needed to obtain the solution, hence state-of-the-art numerical methods are usu- 
ally sought. Recently, the Discrete Geometric Approach (DGA) gained popularity, 
becoming an attractive method to solve BVP arising in various physical theories, see 
for example |4j-|13]. The DGA bears strong similarities to compatible or mimetic dis- 
cretizations [IS] , [IS] , discrete exterior calculus [I?] and finite element exterior calculus 
[18] . [To] , [20] . All these methods present some pedagogical advantages with respect to 
the standard widely used Finite Element Method (FEM). 

First of all, the topological nature of Maxwell's equations and the geometric struc- 
ture behind them allows to reformulate the mathematical description of electromag- 
netism directly in algebraic form. Such a reformulation can be formalized in an elegant 
way by using algebraic topology [5], [6], [8], [9], [16], [18], [I9], [20]. Taking advantage 
of this formalism, as illustrated in Section 3, physical variables are modeled as cochains 
and Maxwell's laws are enforced by means of the coboundary operator. Information 
about the metric and the physical properties of the materials is encoded in the constitu- 
tive relations, that are modeled as discrete counterparts of the Hodge star operator [H] , 
[TT] , [21] . [16] . [20] usually called constitutive matrices [13]. By combining Maxwell's 
with constitutive matrices, an algebraic system of equations is directly obtained, yield- 
ing to a simple, accurate and efficient numerical technique. The difference of the DGA 
with respect to similar methods lies in the computation of the constitutive matrices, 
which in the DGA framework is based on a closed-form geometric construction. For 
a computational domain discretized by using a geometric realization of a polyhedral 
cell complex, one may use the techniques described in [22], [23] and references therein, 
without losing the symmetry, positive-definiteness and consistency of the constitutive 
matrices which guarantee the convergence of the method. Hence, we consider the most 
general situation of dealing with a polyhedral cell complex. 

Our purpose is not to present the widely known DGA or similar discretizations, but 
to use it as a working framework. This choice does not limit the generality of the results. 
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since the standard Finite Element Method (FEM) and the Finite Differences (FD) can 
be easily reinterpreted in the DGA framework as in [8], [10], [11], [2l|, [25], [16], [18], 
[IS], [20]- Consequently our results can be extended, without any modification, to the 
corresponding widely used FEM formulation. 

The paper is focused on a particular application of algebraic topology, namely the 
potential design for the efficient numerical solution of eddy-currents Boundary Value 
Problems (BVP). Electromagnetic potentials are auxiliary quantities frequently used 
to enforce some of the Maxwell's laws implicitly. There are two families of formulations 
for magneto-quasi-static problems, depending on the set of potentials chosen, see for 
example [2], [3], [11]. To better understand the link between (co)homology theory and 
physics, our attention is focused on the /i-oriented geometric formulations, namely the 
T-fl [26], [27], which are based on a magnetic scalar potential Q. Those are much more 
efficient than the complementary family of 6-oriented A and A-x formulations [2H], 
both in terms of memory requirements and simulation time. The main reason is that 
usually /i-oriented formulations require about an order of magnitude less unknowns. 
Nonetheless, when /i-oriented formulations involve electrically conductive regions with 
holes (i.e., the first homology group of some conductor is non-trivial), the design of 
potentials is not straightforward. Cuts are needed to be introduced to make the BVP 
well defined. How to define cuts and devise an efficient and automatic algorithm to 
compute them has been an intellectual challenge for the computational electromag- 
netics community for the last twenty-five years. While a connection of this issue with 
homology theory was quickly recognized by Kotiuga [29] more than twenty years ago, 
heuristic definition of cuts based on intuition are surprisingly still dominant in the 
literature. 

The aim of the paper is to rigorously present a systematic design of the potentials 
employed in /i-oriented formulations by taking advantage of homology and cohomol- 
ogy theories. In particular, at the end of the presentation we are able to formally 
demonstrate that cuts are generators of the first cohomology group over integers of the 
insulating region. The originality of the approach presented in this paper lies also in 
the fact that the design of potentials is tackled directly within a topological setting. In 
fact, thanks to the reformulation of Maxwell's laws by using the coboundary operator, 
homology and cohomology with integer coefficients are employed from the beginning 
for the potential design in place of the standard de Rham cohomology, see for example 
[50] . routinely used in the FEM context, see for example ^SlJ, [29j, [33], [32] • In the 
FEM framework, the so-called non-local basis functions are added to the set of usual 
basis functions to be able to span the de Rham first cohomology group, see for example 
[31], [32], [31], [35], [39], [36], [3Z], [33] • Moreover, employing the DGA, new insights 
into the formulation can be presented by exploiting the dualities arising when, as in the 
DGA framework, two interlocked cell complexes — one dual of the other — are employed. 
For example, the physical interpretation of the non-local basis functions as non-local 
Faraday's equations will become apparent. 
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The second purpose of the paper is to present a survey on definitions of cuts aheady 
presented in the hterature showing their shortcomings. Concrete counter-examples 
show why heuristic definitions should be abandoned and that cohomology is not one 
of the possible options but something which is expressly needed to face this problem. 

The paper is structured as follows. In Section [2} a survey on the relevant topics 
of algebraic topology is provided together with a link to electromagnetic modeling. 
In Section |3| Maxwell's laws casted in algebraic form are recalled. Section |4] shows 
the problem related with expressing Ampere's law with a magnetic scalar potential 
when dealing with conducting regions with holes. Then, to solve this issue, cuts are 
defined as generators of the first cohomology group over integers of the insulating 
region. In Section |5} the T-Q geometric formulation to solve magneto-quasi-static 
BVP is described. Section [6] contains a survey of the definitions of cuts presented 
in the literature together with an illustration of their shortcomings. In Section [7| a 
short discussion on how to compute the cohomology generators is presented. Finally, 
in Section [8| the conclusions are drawn. 

2 Algebraic topology in computational physics 
2.1 Basic concepts in algebraic topology 

In this Section, some basic concepts of algebraic topology are reviewed. Let us first 
introduce the concept of finite regular CW-complex. An n-cell e" is an open subset of 
a Hausdorff space X homeomorphic to the n-dimensional unit ball -B"(0) C M". An 
n-cell e"" is said to be attached to the closed subset K <Z X if there exists a continuous 
map / : -B" (0) — )■ e" such that / maps the open ball -B"(0) homeomorphically onto e" 
and f{dB^{0)) C -ft' in a way that e"ni^ = 0. The map / is referred to as characteristic 
map. The finite CW-complexes are defined in the following way: 

Definition 2.1. Let X denote a Hausdorff space. A closed subset K <Z X is called a 
(finite) CW-complex of dimension N, if there exists an ascending sequence of closed 
subspaces C C . . . C K'^ = K such that the following holds. 

(i) is a finite space. 

(a) For n G {1,...,A^}, the set K"^ is obtained from K"'~^ by attaching a finite 
collection Kn of n-cells. 

In this case, the subset K"^ is called the n-skeleton and the elements of are 
called the vertices of K. An A^-dimensional CW-complex is called regular if for each 
cell e", where n G {1, . . . , A^}, there exists a characteristic map / : -B"(0) — )■ e" which 
is a homeomorphism on i?"(0). In this case, we say that the m-cell e"* is a face of an 
n-cell e", if the inclusion C e" holds, see [ID]. Moreover, since our aim is to model 
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physical objects, we restrict to the case of regular CW-complexes embedded in M^. 
Finally, by a polyhedral mesh we mean a cellular decomposition of the considered space, 
which is a regular CW-complex such that each cell of the complex is a polyhedron. In 
the paper we are going to use the terms mesh and regular CW-complex interchangeably. 

Let /C be a collection of cells of the regular CW-complex K. Let k : /Cj_i x K-i ^ 
{ — 1,0,1} for z G Z be the so-called incidence index which assigns to a pair of cells 
their incidence number (for further details consult |1Q]). Let G denote the module of 
integers (Z), real (M) or complex (C) numbers. The group of formal sums YleeK "-^e^' 
where E G for every e G /Cj, is the group of i- chains of the complex /C and is denoted 
by Ci{)C,G). For a chain c = J2eejc,^ee the support |c| of c consist of all elements 
e G /Cj such that 7^ 0. For two chains c = ^eeA:^ and d = J2eejc, their scalar 
product is (c, c?) = Xleex:, '^e/3e- The group of cochains C*(/C,G) is formally defined 
as the group of maps from elements of Cj(/C, G) to G with coordinatewise addition. 
However, it is possible and convenient for the computations to represent a cochain as 
a chain. Namely, to determine the value of a map c* : Ci(/C, G) — )■ G on any z-chain, it 
suffices to know the value of c* on every e G /Cj. In this way, it is possible to associate 
to the cochain c* a chain c such that for any other chain d G Ci{IC,G) the value of 
cochain c* on chain d is equal to (c, d). For a cochain c* its support \c*\ consists of all 
the cells whose value of c* is nonzero. 

In this paper, two kind of cochains are considered. The first are the integer-valued 
cochains — for example, the representatives of the first cohomology group generators 
over integers. The second are the complex-valued cochains, which model physical 



variables in the proposed application as discussed in Section 2.3 

Let us define the boundary map di : Ci{)C, G) — )■ C'j_i(/C, G). For an element e G /Cj 
we define 

die= ^ K(/,e)/ 

and extend it linearly to the map from Gi{IC,G) to C'i_i(/C, G). The coboundary map 
6' : G%JC, G) G'+\JC, G) is defined for e G Id by 

and extended linearly to the map from G*(/C,G) to G*"'"^(/C, G). It is standard that 
di-idi = 5*(5*~^ = for every i G Z, see [H]. Moreover, the coboundary map is dual to 
boundary map in homology. In fact, it can be equivalently defined with the equality 
{6c*, d) = {c\ dd) for every c* G G'-\K, G) and for every G G,(/C, G), see [H]. 

The boundary operator gives rise to a classification of chains. The group of i- 
cycles is Zi{fC,G) = {c G Gj(/C,G)|9c = 0}. The group of i-houndaries is Bi{]C,G) = 
{c G Gi(/C,G)| there exist d G Ci+i{}C,G)\dd = c}. Intuitively, a cycle is a chain 
whose boundary vanishes while a boundary is a cycle which can be obtained as the 
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boundary of some higher dimensional chain. The ith homology group is the quotient 
Hi{]C,G) = Zi{IC,G)/Bi{}C,G). The cycles that are not boundaries are nonzero in 
Hi{)C,G). The cycles that differ a by a boundary are in the same equivalence class. 
Given a chain c, by [c] we denote its homology class, i.e. the class containing all the 
cycles homologous to c. By generators of the ith homology group we mean a minimal 
set of classes which generates Hi{}C,G). In the following, for the sake of brevity, by 
generators we also mean the cycles being representatives of the considered classes that 
generate Hi{]C,G). 

Dually, with the coboundary operator the cochains may be classified. The group 
of i-cocycles is Z^{}C,G) = {c G C*(/C,G)|(5c = 0}. The group of i-coboundaries is 
B\1C,G) = {c G C\1C,G)\ there exist d G C'-\lC,G)\5d = c}. The ith cohomology 
group is the quotient H\}C,G) = Z\)C,G)/B\}C,G). By generators of the ith co- 
homology group we mean a minimal set of classes which generates H^{]C,G). Also 
in this case, for the sake of brevity, by generators we also mean the cocycles being 
representatives of the considered classes that generate W{]C,G). 

In the following, we will use also the standard concept of the so-called relative 
(co)homology. In relative (co)homology, some parts of the complex may be consid- 
ered irrelevant. Let K, be the considered regular CW-complex and 5 C /C be a closed 
sub-complex of JC. The concept of relative homology bases on the definition of relative 
chains Ci{)C,S,G) = Ci{}C,G)/Ci{S,G). The definition of relative cycles Zi{]C,S,G), 
relative boundaries Bi{fC,S ,G) and relative homology group Hi{}C,S,G) remain un- 
changed with respect to the absolute version once relative chains are used. Exactly the 
same approach is also employed in defining relative cohomology. 



In theorem 2.3, it will be recalled that there is no torsion |H] in the homology 
and cohomology groups dealing with regular CW-complexes embedded in M^. A direct 
consequence of the Universal Coefficient Theorem for cohomology, see [H], is that in 
the considered torsion-free case the generators of the cohomology group over integers 
and the generators of the cohomology group over complex numbers are in a bijective 
correspondence (for further details see |32])- Therefore, all the computations are rigor- 
ously performed by using integer arithmetic and the obtained cohomology generators 
are valid cohomology generators also in the case of complex coefficients. The theory of 
(co)homology computations for regular CW-complexes can be found in |10]. The impor- 
tant point is that for (co)homology computations only the incidence indexes k between 
cells of K are needed. This fact provides an easy way of representing CW-complexes 
for the (co)homology computations with a computer by using a pointer-based data 
structure. Once the theory is provided, any of the existing libraries like [35] or |16] 
used to compute homology of, for instance, cubical sets can be adopted for cohomology 
computation of an arbitrary regular CW-complex. 

Let us now introduce the concept of exact sequences. Let Ai, . . . , A^+i be abelian 
groups and let : — )■ Aj+i for i G {!,... ,m} be homomorphisms between them. 
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The sequence 



A ^ A ^ ""'-I, A ^ A 
/il > /i2 7- . . . > Am r /im+1 



is called an exact sequence if Im(aj) = Ker(ai+i) for every i G {1, ... ,m — 1}. For 
m = 4 and = ^4 = the exact sequence 

is referred to as short exact sequence. The so-called exact sequence of the reduced 
homology provides us a tool to relate the homology group of the space X, its subspace 
A and the relative homology of the pair {X, A). 

Theorem 2.2 (see Theorem 2.16, [H]). If X is a regular CW-complex and A G X is 
a sub-complex of X , then there is a long exact sequence 

... 4 Hr,{A) ^ Hn{X) ^ Hn{X,A) 4 ^ . . . . 

The map d : Hn{X,A) — )■ maps a class [a] G Hn{X,A) to a class [da] G 

Hn-i{A). It is straightforward that, when Hn{X) is trivial, from the exactness of the 
sequence, d : Hn{X,A) — )■ if„_i(y4) is an isomorphism. 

To state, in Section [4], the dualities between first homology and first cohomology 
groups of subsets of the following theorems are required: 

Theorem 2.3 ([H], Proposition A.4, Corollary 3.44). // X C zs a fimte CW- 
complex, then Hi{X, Z) is for i > n and torsion free for i = n — 1 and i = n — 2. 

For n = 3 the Theorem 2.3| states that the first and the second homology group of 
a finite CW-complexes embeddable in are torsion free. 

For CW-complexes the following — stronger from standard — version of excision the- 
orem holds: 

Theorem 2.4 (Corollary 2.24, pfT]). // the CW-complex X is the union of sub- 
complexes A and B, then the inclusion {B, An B) "-^ {X, A) induces an isomorphism 
Hn{B, A n 5) ^ Hn{X, A) for all n. 



2.2 Dual chain complex 

A polyhedral mesh is used to model the domain of interest of the electromagnetic 
problem. Let us fix the polyhedral mesh /C = {Kn}n&i- Let us now define the dual 
mesh B. The construction is a straightforward extension of the construction of the 
dual mesh for a simplicial complex explained in Figure [l| see [H] . 

Let no = dim /C be the dimension of the complex /C. For every cell c G Ki for 
i G {0, . . . ,no} by c G B^o-i let us denote the corresponding element in dual mesh B. 
For every c G KnQ, the corresponding element c G is simply the barycenter of c. 
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a) b) c) d) 



Figure 1: a) A cell f of a simplicial complex and its one-to-one node of the barycentric 
complex ng; b) A face / of the simplicial complex and its one-to-one edge of the 
barycentric complex eg; c) An edge e of the simplicial complex and its one-to-one face 
of the barycentric complex /g; d) One-to-one correspondence between a node n of the 
simplicial complex and the volume of the barycentric complex wg. 

The remaining cells of B are defined recursively in the following way. For c & Ki for 
i G {0, . . . , no — 1}, let {ci, . . . , c„} = \5c\ and let -B(c) denotes the barycenter of c. 
Then c = IJiLi Uzeciit-^' -^('^)]}) where [x, y] denotes the line segment joining x and y. 
In the paper, we denote by d and 5 the boundary and coboundary operator in the dual 
complex B. To let the chain complex of B be dual — in purely algebraic sense — to the 
chain complex of /C, the boundary operator is defined as follows: 

(5c, d)K. = {dd, c)b Vc G Ci(/C), G Ci_i(/C) for i G N. 

We would like to point out that the presented construction is valid only in the case of 
manifolds without boundary. In case of the presence of a boundary, first the complex 
dual to the boundary is constructed. Then, the two dual complexes, namely the one 
described above and the complex dual to the boundary, are merged in an obvious way. 
Nonetheless, the construction of the dual complex is fundamental only to develop the 
formulation while, for the computations, the explicit construction of the dual complex 
can be avoided. The historical context of the idea of barycentric dual complex is 
mentioned in Section [4.3. II 

2.3 Physical variables as complex- valued cochains 

Let /C be a homologically trivial polyhedral cell complex in M"^. Let /Cc be a closed 
sub-complex of /C which models the conducting region. The K, \ int ICc is an insulating 
region the closure of which is meshed by /Cq, which is a closed sub-complex of /C. The 
interface of the conducting and insulating region is meshed by /Cc n /C^. Moreover, it 
is assumed that /Cc and )Ca are non-empty. The dual sub-complexes are denoted by Ba 
and Be, respectively. 

Since our aim is to solve the eddy-current problem in the frequency domain, the 
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physical variables are modeled in this paper as complex-valued cocliainq^ According 
to Tonti's classification of physical variables [6], [9], [13], there is a unique association 
between a physical variable, such as electric current or magnetic flux pQ, [B], [H], and 
an oriented geometric element of the two cell complexes /C and B. The cochain values, 
usually called Degrees of Freedom (DoFs) in computational physics (see for example 
[H]), have a direct physical interpretation: By using the so-called de Rham map [50] , 
they are defined as integrals of the electromagnetic differential forms over the elements 
of the complexrj 




Figure 2: Association of the degrees of freedom to the oriented geometrical entities. 

The focus of this paper is on /i-oriented formulations [2S], [2Z|, so the following 
physical variables and related association with geometrical elements of /C or i3 are 
considered: 

• (I, /) is the electric current associated to the face f E JC, see Fig. |2^. (I, /) = 
over the faces / G /Ca (with this definition, the current associated with the faces 
/ e /Ca n /Cc is set to zero, since there is the need of a boundary condition that 
prevents the current to flow thought the boundary of the conductive region); 

• (F, e) is the magneto-motive force (m.m.f.) associated to the edge e G /C, see Fig. 

• /b) is the magnetic flux associated to the dual face fs G B, see Fig. [ij^; 

• (U, eg) is the electro-motive force (e.m.f.) associated to the dual edge G B, 
see Fig. |2]i. 

frequency domain eddy-current problem implies that all physical variables exhibit a time vari- 
ation as isofrequential sinusoids. By using the standard symbolic method, see |49j . constant complex 
numbers called phasors are used to represent the sinusoids. If the eddy-current problem has to be 
solved in time domain, the reals should be used in place of the complex numbers through the paper 
without any further modification. 

^For example, the magneto-motive force (m.m.f.) DoF relative to the 1-dimensional cell e is the 
integral of the differential 1-form magnetic field over e. 
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All these complex values, one for each 1- or 2-dimensional cell in the corresponding 
complex, are the coefficients of the corresponding complex-valued cochains denoted in 
boldface type. For a fixed chain and cochain basis, each cochain can be represented as 
a vector which is used in the computations. In the following, since no confusion can 
arise, the notations for (co)chains and vectors representing them will be the same. 

3 Maxwell's equations in algebraic form and poten- 
tials analysis 

In this Section, the algebraic Maxwell's laws |5], [6], [7] are reviewed. 

The discrete current continuity law enforces the dot product of the currents asso- 
ciated with faces belonging to the boundary of a volume v G /Cc, with dim f = 3, to 
be zero 

{I,dv) = {6I,v) = 0, Vt; G ICc. (1) 

Focusing on the generic face /, the discrete Ampere's law enforces the dot product 
of the m.m.f. F on the boundary of the face / to match the current associated with /, 

(F,9/) = (5F,/) = (I,/), V/G/C. (2) 

Since (I, /) = 0, V/ G ICa, F is a 1-cocycle in /C^ (however, it is not a cocycle in /C). 

The discrete magnetic Gauss's law enforces the dot product of the magnetic fluxes 
associated with the dual faces belonging to the boundary of a dual volume to be 
zero 

{^,dvB) = {6^,VB) = 0, yvsEB. (3) 

Focusing on a dual face fs, the discrete Faraday's law enforces the dot product 
of the e.m.f. U on the boundary of /g to match the opposite of the variation of the 
magnetic flux through the face /g. Considering problems in frequency domain, this 
translates in 

(U, dfs) = (5U, h) = fs), e B, (4) 

where u is the angular frequency of the sinusoids equal to 2tt times the considered 
frequency. 

The expressions just discussed of the four algebraic laws are called 'local'. There 
exist also the so-called 'non-local' versions of each of them, which are obtained by 
considering the balance not on exactly one geometrical entity but on a chain. It is 
straightforward to see that if the region is homologically trivial each non-local law can 
be obtained by considering a linear combination of local laws. Therefore, in this case, 
the non local laws do not bring any new information. As will be discussed in the next 
Sections, this does not hold for homologically non-trivial regions. 
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3.1 A preliminary definition of potentials 

In this Section a preliminary definition of potentials employed in the formulation is 
presented and analyzed by means of algebraic topology. For this preliminary definition 
let us assume that the considered complex /Cc is homologically trivial. In particular, 
we use the fact that, when the conductors and the whole domain are homologically 
trivial cell complexes then, from a standard result on exact sequences, one has that 
the homology of the insulating domain is also trivial. How to generalize the definition 



of potentials in case of homologically non-trivial regions is the subject of Section 4.1 

Let us first analyze the potentials employed in the insulating region ICa- It has 
been already indicated that F is a 1-cocycle in /C^, hence (5F, /) = holds V/ G ICa- 
Thus, since ICa is homologically trivial, a magnetic scalar potential 0-cochain can be 
introduced in the insulating region such that 

(5n,e) = (F,e),VeG/C,. (5) 

Let us analyze now the potentials employed in the conducting region /Cc. From 0, 
we know that I is a 2-cocycle in /Cc, hence (51, f) = holds Vt; G /Cc- Thus, an electric 
vector potential T 1-cochain can be introduced in the conducting region such that 

(T,9/) = (5T,/) = (I,/),V/G/Cc. (6) 

Thanks to (|2]), the following holds 

(5T,/) = (I,/) = (5F,/),V/G/Cc. (7) 

Since F and T are 1-cochains such that 5F = 5T, they differ by a 0-coboundary of a 
0-cochain O, i.e. F = T + 50. Since it is required that the magnetic scalar potential 
is continuous inside /C, we can extend the support of Vt also inside /Cc in such a way 
that (O, n) = (O, n) for every node n G /Cc. We want to remark, that those extensions 
are valid also in case of homologically non-trivial complexes /Cc and /C^ and they are 
used further in the paper. For brevity, let us define also the cochain T as cochain in 
/C. To this aim, we assume that (T, e) = for every edge e G ICa- 

To sum up, by using the potentials T and fl, Ampere's law ^ and current con- 
tinuity law ([!]) can be enforced implicitly by considering the following expression for 
F 

F = 6n + T. (8) 

Then, it is easy to show that employing this definition of potentials Ampere's law 
holds in ICa (in fact, {6FJ) = {65nj) = (0,/) = = (I,/), V/ G JCa, since the 
current is zero in ICa) and the current continuity law holds in all /C (in fact, SI = SST = 



0). The remaining laws will be enforced by a system of equations in Section 5.2 
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4 Potentials design 



4.1 Non-local Ampere's law in homologically non-trivial do- 
mains 

Let us now remove the hypothesis of /Cc being homologically trivial. Ampere's law can 
be written on a 1-cycle c G Zi(/Ca) and a 2-chain s G C2(/C) such that ds = c: 

(F,c) = (I,s). (9) 

Equation ^ is an example of a non-local equation, since the algebraic constraint is not 
enforced on a geometric element, like (|2|, but involves geometric elements belonging 
to a wider collection, in this case the support of s and its boundary c. In ^ we have 
not specified which 2-chain s has to be used for taking the dot product at the right- 
hand side (in fact, s has been determined only up to its boundary ds = c). To solve 
this issue, we show with next Theorem that the value (I, s) depends only on ds and, 
therefore, ^ is well defined. 

Theorem 4.1. Let si and S2 he two 2-chains such that dsi = ds2 = c, where c G 
Zi(/Ca). r/ien (I, si) = (I, S2). 

Proof. From the assumptions, si — S2 G Z2{)C). Since /C is homologically trivial, there 
exists b G C^^K,) such that db = si — S2- Consequently, si = S2 + db holds. Then, 
(I, si) = (I, S2) + (I, db) = (I, S2) + {51, b) = (I, S2), since, due to Q, 51 = 0. □ 



Due to Theorem 4.1[ one can state the following definition: 

Definition 4.2. The value Ic = for an arbitrary 1-cycle c G Zi{ICa) and a 

2-chain s such that c = ds, is called current linked by the 1-cycle c. 

Let us now show that the current linked by a 1-cycle c G Zi{ICa) is the same for all 
the cycles in the homology class of c. 

Theorem 4.3. The linked current Ic = (I, s) depends only on the ifi(/Ca) class of the 
1-cycle c = ds, where c G Zi(/Ca). Thanks to non-local Ampere's law, the same holds 
for the dot product (F, c). 

Proof. Let us take two cycles Ci, C2 G /Ca in the same homology class ifi(/Ca). It means 
that ci = C2 + ds holds, where s G C2{ICa). Then, (F,ci) = (F,C2) + {F,ds) = 
(F, C2) + {SF, s) = (F, C2) + (I, s). From the definition of I we have that (I, /) = for 
every / G /Ca, therefore (I, s) = 0, what completes the proof. □ 
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The idea of the proof of Theorem 4^ is presented in Fig. |3] for the complement 
of a sohd double torus with respect to a cube which contains it (not represented in 
the picture for the sake of clarity). The supports of two cycles Ci and C2 in the same 
homology class are depicted in Fig. |3^. The supports of two 2-chains Si and S2 used 
to evaluate the linked currents 1^ = (I, si) and Ic^ = (I, are depicted in Fig. |3}d. 
Finally, in Fig. ^ it is possible to see the support of a 2-chain s such that ci = C2 + ds, 
which has been used in the proof. 




Figure 3: a) Two cycles, Ci and C2, in the same HiiJCa) class, b) Example of two 
possible 2-chains used to evaluate the currents linked to ci and C2. c) The 2-chain s 
whose boundary 1-chain is ci — C2. 

Using the scalar potential in /Cq, as the definition ([s]), yields to an inconsistency, 
since Ampere's law may be violated on some 1-cycle c = ds E Zi{)Ca)- In fact, let us 
consider a 1-cycle c G ifi(/Ca). Using the double torus example previously introduced, 
this cycle may be for example the cycle ci or C2 in Fig. [3^. Now, due to (|8]), we have 

(F, c) = (F, ds) = (T + 6n, ds) = (T, ds) + {6n, ds) = {66n, s) = (0, s) = ^ (I, s), 

since (T, e) = for every e G )Ca- The last inequality 7^ (I, s) follows from the fact 
that the current flowing through s is non-zero in general, since the support of s has to 
intersect /Cc- According to next Theorem, the presented problem may occur only with 
cycles which are non-trivial in HiiJCa). 

Theorem 4.4. Ampere's law is satisfied for every 1-boundary b in JCa- The cycles that 
produce an inconsistency in Ampere's law — because may link a non-zero current — are 
the cycles which are non-trivial in the 1-st homology group Hi{lCa)- 

Proof. In the first case, the cycle c G Bi{)Ca) is bounding, which means that there exists 
s G C2(/Ca) such that ds = c. Since s does not intersect /Cc, the dot product of the 
current I and s is zero and the dot product of the m.m.f. F and c is also zero. 

In the second case, c being non-zero in the first homology group Hi{)Ca), such a 
chain s G C2{JCa) does not exist. But, since the whole complex /C is homologically 
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trivial, there exist a 2-chain s' G C2{1C) such that ds' = c. Consequently one needs to 
have \s'\ fl/Cc 7^ 0- Thus, the support of s' has to extend in the current- carrying region 
/Cc and, consequently, the current trough s' (and linked by c) is non-zero in general. 
This causes in general an inconsistency in Ampere's law for the cycles c non-trivial in 

Hl{lCa). □ 




Figure 4: a) Two non-trivial elements of i7i(/Ca) called ci and C2. b) Let us consider 
2-chains in JC whose boundaries are ci and C2- c) Two independent currents, ii and i2, 
evaluated through 2-chains si and S2 or through 2-chains cxi and a2- d) The 1-cycles 
dai and da2 are in the same homology HiiJCa) class as the corresponding 1-cycle Ci 
and C2- This is formalized in the following of the paper. 

The idea of last proof can be presented by using again the double torus example. 
Let us introduce two non-trivial elements of Hi{lCa) called Ci and C2 which are the 
representatives of generators of the first homology group HiifCg) and whose support 
is depicted in Fig. |4^. The support of two 2-chains Si, i G {1,2}, such that dsi = Ci, 
are represented in Fig. |4]d. It is easy to see that the supports of these 2-chains have to 
intersect /Cc- 
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4.2 Independent currents 

After showing the inconsistency arising in general from the definition of potentials 



presented in Section 3.1, some modifications in this definition are needed for Ampere's 
law to hold implicitly for every 1-cycle c G Zi{ICa) and a 2-chain s such that c = ds. 
We define in this Section a key ingredient for the new definition of potentials, which is 
called set of independent currents. 

It has been already pointed out that the electric current is nonzero inside the 
conducting region ICc only. Consequently, for every chain s G C2(/C) for which ds = 
c G Zi{ICa), the current can be non-zero only on |s| fl /Cc. Since ICc is a sub-complex of 
/C, the restriction of the 2-chain s to /Cc, denoted as a = s\ic^, is an element of C2{JCc)- 
The supports of such restrictions of the 2-chains si and S2 in the double torus example 
are the cxi and (T2 shown in Fig. Since da G Ci(9/Cc), we have that a G Z2{}Cc, 9/Cc). 
Consequently a can be generated from H2{}Cc, <9/Cc) basis by adding the boundary of a 
suitable 3-chain d G C^llCc)- The current through any chain non-zero in H2{K,ci dJCc) is 
determined by the current through the H2{)Cc, dlCc) basis elements (for trivial 2-chains 
and relative 2-chains the current is zero as a consequence of local version of current 
continuity law ([T])). 

For the whole Section, let us fix the set of relative cycles ai, . . . , cr„ G Z2(JCc, d)Cc) 
representing the homology group H2{}Cc,d}Cc) generators. From the Theorem 2.2 one 
may write the following exact sequence of the pair (/C, /Ca): 

... 4 H2i}Ca) ^ H2{}C) ^ H2{}C, Ka) A H^{Ka) h.... 

The assumption that the mesh K, is acyclic provides H2{1C) = 0. Consequently, d : 
H2{}C,)Ca) HiiJCa) is an isomorphism defined in the following way d : H2{}C,}Ca) 9 
[a] ^ [da] G H.iKa). 



Let us use Theorem |2.4| for the sub-complex A equal to /Cq and sub-complex B 
equal to /Cc and X equal to /C. Since ICaP^lCc = dJCc, the following inclusion map 
(/Cc, 5/Cc) 7- (/C, /Ca) induces the isomorphism if2(/Cc, /Cc fl /Ca) = H2{JCc,d}Cc) 
i72(/C,/Ca). Consequently, we have the following sequence of isomorphisms 

/f,(/C„9/Cc) //2(/C,/Ca) A i/i(/Ca), 

where on the left-hand side there is the group generated by the if2(/Cc, dlCc) generators 
and on the right-hand side there is the group generated by the classes of cycles on 
which the Ampere's law has to be enforced. The image of the generators of if2(/Cc, OK-c) 
through the above isomorphism is referred to as independent cycles. 
What we showed in this Section motivates the following definition: 

Definition 4.5. The complex numbers being the dot-products of the current I with the 
representatives of a basis {crj} of H2{ICc,dlCc) are called independent currents {ij} in 

/Cc 

ij = {I,(Tj),J G {l,...,/3i(/Ca)}. 
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We want to point out that the {aj} are integer H2{JCc,d}Cc) homology generators. 
When evaluating the dot-product, we trivially interpret them as complex homology 
group generators. 

Since the isomorphism H2{}Cc, dlCc) -f^2(^, ^a) is induced by the inclusion map, 
ai, . . . , (T„ form also a set of cycles representing the basis of H2{}C, K-a)- When passing 
by the isomorphism induced by the boundary map d : H2{}C, )Ca) — ifi(/Ca), one gets 
that dai, . . . , 9cr„ is a set of cycles representing a basis in Hi{lCa). 



4.3 Definition of cuts 

First of all, let us note that it suffices to enforce Ampere's law on the cycles 
[5(Ti], . . . , [5(T„]. Then, for every other cycle [c] G HiiJCa) such that [c] = [^"^^ AjScXj], 
the current linked by c is equal to Y^^=i^i^i^i) — Sr=i -^^^i' which follows from the 
fact that dot product of the m.m.f. on boundaries is zero. 

Now, taking into account the arguments in the last Section, we would like to modify 
the definition ([s]) of F in K,a in such a way that Ampere's law is satisfied for all cycles 
c G ZiiJCa)- Since F is a 1-cocycle in /Ca, we are going to construct a family of 1- 
cocycles {c*}^^^ in Ka over Z which, after being multiplied by the independent currents 
are added to T. In particular, the family of 1-cocycles {c*}"^^ should verify 
(c*, daj) = 6ij for every i, j G {1, . . . , n}. We are going to show that, for this purpose, 
the representatives of a basis of the 1-st cohomology group H^{lCa) dual to the HiiJCa) 
basis [c?o"i], . . . , [ScXn] are needed. To prove the existence of the dual basis, let us recall 
the Universal Coefficient Theorem for cohomology. 

Theorem 4.6 ([S], Theorem 3.2). // a complex K,a has (integer) homology groups 
Hn{lCa), then the cohomology groups H'^{)Ca,G) are determined by splitting exact se- 
quences 

^ Ext{Hn-iilCa), G) H'^iJCa, G) A HomiH^JCa), G) ^ 0. 



In this paper, there is no need to go into the definition of the Ext functor. The key 
property is that Ext{Q, G) = if Q is a free group. For further details and proof of 
this property consult [H]. 

For a class [d\ G H"-{}Ca,G), since is a cocycle, one has = {6d, z) = {d,dz) for 
every z G Z„_|_i(/Ca). From the above equality, it follows that d\B„(iCa) = 0- Let us define 
the restriction do = d\z„(Ka)- Since (io|B„(K:a) = 0, then d^ G Hom{Hn{}Ca),G). This 
shows the correctness of the definition of the map h{[d]) = do E Hom{Hn{}Ca), G) used 
in the exact sequence in Theorem 4.6 Finally, we need to show that, for d, d' G [d], one 
has h{d) = h{d'). Since d,d' G [rf], there exists e G C^{K) such that d = d' + 6e. Let 
us take a cycle / G Ci(/C), then we have df = 0. Consequently {d, f) = {d' + 6e, f) = 
{d', f) + {Se, f) = {d', /). Therefore, the value of the map h does not depend on the 
representatives of the cocylce and h is well defined. 
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In our case the group G is the group of integers and the Universal Coefficient 
Theorem for cohomology is used for n = 1. In this case, the exact sequence from 
Theorem 14.61 has the form 

^ Ext{Ho{JCa),Z) H\JCa,Z) ^ Hom{Hi{ICa),Z) 0. 

For the complex ICa we have that ifo(A^a) = for some p G Z, p > 0. This 
provides that Hq^JCo) is a free group. From the cited property of the Ext functor, it 
follows that Ext{HQ{ICa),'^) = 0. From the exactness of the sequence, one has that 
h : i7^(/Ca,Z) — -f^om(/7i(/Ca), Z) is an isomorphism. 



Due to Theorem 2.3 the homology group Hi{ICa) is torsion free. This provides, from 
Theorem 3.61 in [52], that it is isomorphic to the direct sum of dim(ifi(/Ca)) = n copies 
of Z. From the set of cycles dai, . . . ,dan forming a HiiJCg) basis, a set of functions 
Ci, 2 G {1, . . . ,n} such that C,i{\daj\) = 6ij form a basis of ifom(ifi(/Ca), Z). From the 
description of the isomorphism h : H^iJCa, Z) — > Hom{Hi{}Ca), Z), it is straightforward 
that h-^iCi) is a cochain c* being an element of the H^{)Ca, ^) basis we are looking for. 

In the presented reasoning we have started from the set of independent currents 
and end up to the cohomology basis. The Reader should be aware that exactly the 
same reasoning can be made the other way around. Namely, if one starts from the 
H^iJCa) basis, then is able to find the corresponding H2{JC,}Ca) basis, which directly 
corresponds to some basis H2{lCc,dlCc) yielding a set of independent currents. 

The following Theorem follows easily form what has been already said. 

Theorem 4.7. Let {c*}fj^^/^°-* be the cocycles representing the H^{ICa) basis. Let 
{dai}^!}^''^ be the cycles representing the dual Hi{JCa) basis. Once we redefine the 
m.m.f. as F = 5Q, + T + X^fii'^"'* ^^en the current linked by the cycles in the ho- 
mology class of dai is equal to ii. Ampere's law holds for every 1-cycle c G Zi{}Ca). 
Hence, the potentials are now consistently designed. 

Further in this paper we assume that the cocycles {c*}f^^/^'"'* representing the H^{K,a) 
basis are defined in the whole complex /C. Therefore, we assume that (c*, e) = for 
every i G {1, . . . , /3i(/Ca)} and for every e G /Cc \ ^a- 

Thus we are now able to give the definition of cuts: 

Definition 4.8. The cuts {c^}^^}^ " are defined as representatives of the first cohomol- 
ogy group generators over integers of the insulating region /Cq. 

To illustrate the presented idea let us consider the two fixed generators Ci and 
C2 for Hi{K,a) relative to the previous example, see Fig. |5^. In the same Figure, a 
concrete example of a representative of a H^(JCa) cohomology generator — dual to the 
HiiJCa) generator ci — is shown. The edges in the picture are the ones that constitute 
the support of c^. The integers associated to these edges are given in such a way that 
{c^,Ci) = 1 and {0^,02) = 0. 
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Figure 5: a) The two generators Ci and and the support of the 1-cochain is shown 
for the double torus example, b) The dual faces /g, dual to edges in the support of c^, 
from a 2-chain Si on the barycentric complex. 



4.3.1 Interpretation of cuts on the dual complex 

Some results presented in this Section are part of an old-fashioned proof of Poincare- 
Lefschetz duality for manifolds with boundary. For this paper, the following restricted 
version of the duality is used: 

Theorem 4.9 (Poincare-Lefschetz duality). H^{Kia) = H2{Ba,dBa). 

The modern proof of this famous Theorem bases on the idea of cup product, as 
for instance in [H]. However, the original proof proposed by Poincare himselQ is 
based on the concept of dual cell structure, which has been described in this paper in 



Section 2.2 For the classical proof of Poincare-Lefschetz duality, one may consult |18] 
or [51]. In this proof the dualization operator D is defined on the complex /Cq in the 
way that, for a 1-cell c G /Ca, the corresponding dual 2-cell Dc G Ba is assigned. The 
presented map turns out to induce the isomorphism in the Poincare-Lefschetz duality 
(for further details consult |18]). 

Consequently, form the Poincare-Lefschetz duality, once the 1-cocycles that repre- 
sent a basis {c^, . . . , c"} of H^iJCa) are provided, it is clear that the set of dual 2-cycles 
di, . . . ,dn defined in the following way: 

di=J2 (^tS, where at = {c\D-^S) 

are the relative cycles that represent a basis of H2{Ba,dBa)- These cycles are denoted 
as {Sjjfi^^/^'"'*. The visualization of the presented duality for the proposed example can 

•^Which turned out not to be complete, but was corrected later on by the Poincare's successors. 



20 



be seen in Fig. [5j On the left, the cohomology generator c* is represented, while, on 
the right, the representative Sj of the generator of H2{Ba,dBa) is depicted. 

5 T-Q magneto-quasi-static formulation 

After the potential design, in this Section we analyze how to solve the magneto-quasi- 
static BVP. 

5.1 Constitutive matrices 

The discrete counterparts of the constitutive laws links fc-cochains in /C with (3 — k)- 
co chains in B: 



* = /xF (a), U 



Q I|yc. (b)- (10) 



The constitutive matrices provide a relation between cochains on the primal and 



cochains on the dual complex (see Section 2.2). The square matrix /j, is called per- 
meance matrix and is the approximate discrete counterpart of the constitutive relation 
B = yU H at continuous level, being the permeability assumed element-wise constant 
and H and B are the magnetic field and the magnetic flux density vector fields, re- 
spectively. The square matrix q is called resistivity matrix and is the approximate 
discrete counterpart of the constitutive relation E = ^ J at continuous level, g being 
the resistivity assumed element-wise a constant and E and J are the electric field and 
the current density vector fields, respectively, g is defined to be zero for geometric 
elements in )Ca \ /Cc- 

Describe in detail how to construct the constitutive matrices g and /j, goes beyond 
the purpose of this paper. Methods valid for a general polyhedral mesh are described 
for example in [22] , [23] and references therein. 



5.2 Algebraic equations 



In this Section, the constitutive matrices described in the Section 5.1 are combined 



with the local algebraic laws presented in Section [3] to obtain an algebraic system of 
equations. 

Up to now, we know that Ampere's law holds in JCa and current continuity law 
holds in /C. Hence, we have to enforce the other laws by means of a linear system of 
equations. 

To do this, let us start from the magnetic Gauss's law (jsj) {6^,vts) = and let us 
use the constitutive relation (lOi) $ = //F. Consequently, we get {6 i-i F, vq) = 0. The 
definition of potentials F = 5Q, + T + X]f=f'"'' ^i^"* ^^om Theorem 



4.7 



is substituted in 
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the last equation. In this way, the final equation is obtained: 



0. 



(11) 



Now, Faraday's law has still to be enforced in the conducting region (due to the 
boundary condition (T, e) = on edges e G 9/Cc, the e.m.f. is not needed in the 
insulating region by using this formulation, see for example [21 p. 1030] for a more 
detailed explanation). To do this, let us start from the local Faraday's law (5U + 
iu^,fs) = 0, V/g G Be- Now, let us substitute the constitutive relations (10) in the 
above equation. In this way, we obtain {6q1 + iuf-iF, f^) = 0, V/g G Be- Let us now 
use the local Ampere's law ^ in JCc by substituting 5F = I 

C6Q6F + tiOfiFJs) = 0,yfB^B,. 

For the sake of brevity, let us define a matrix R = 6q6 + icon. Hence, we can write 



Since R^il = iuf-iSfl, we obtain: 



ti6n + RT + ^ RcHj, h) = 0, V/b G Be. 



(12) 



The unknowns are the coefficients of the cochain f2 associated with each node and the 
coefficients of T associated with edges belonging to /Cc \ JCa (in fact, (T,e) = 0, Ve G 



)Ca). We have written the equations (11) and (12) corresponding to these unknowns. 
But we have also the independent currents ij as additional unknowns. Which are the 
needed additional equations and where do they come from? 



5.3 Non-local Faraday's equations and the final linear system 
of equations 

The dot product of the e.m.f. (U, &) with every bounding 1-cycle b G Ci{Bc) is easily 
determined by using a non-local Faraday's law 

(U, b) = {-ico^, s), be Bi{Bc) and b = ds, 



which is a linear combination of local Faraday's laws already e nforc ed by ( 12 ). Similarly 
to what developed about the independent currents in section 4.2 , the linked fiux <l>c = 
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($,s), linked by the cycle b, does not depend on the 2-chain s. This is because the 
local Gauss's magnetic law ^ hold thanks to (11). Hence, such non-local equations 



written on boundaries do not bring any new constraint. 

On the opposite, {(J,h) over a 1-cycle h G Ci{Bc) nonzero in HiiBc) cannot be 
determined by using only cochains in Be- This is because (U, h) depends on cochains 
in Be and Ba through the non-local Faraday's law. Namely (U, h) has to match the 
magnetic flux variation (—20;$, s) through a 2-chain s such that c = ds. The key point 
is that the support of s extends also in the sub-complex Ba- 

We need to show that the 2-chains used to take dot product with the magnetic flux, 
whose boundary are the Hi{Bc) generators, are generators for H2{B,Bc). Similarly to 
what done with the currents, we need just the H2{Ba, dBa) generators and the reasoning 
presented in this Section is analogous to the one presented in Section |4.3[ Let us take 
a fixed set of 1-cocycles {c^, . . . , c"}, d G C^(/Ca) for i G {1, . . . , n}, representing the 
basis of H^iJCa)- The set of cycles {di, . . . , d^}, di G C2(i3a, dBa) for i G {1, . . . , n}. 



defined in the Section 4.3.1, represents a basis of H2{Ba, dBa)- Let us write the long 



exact sequence of the pair {B, Be)- 

... 4 ^ H„iB) ^ H^B, B,) 4 Hn-iiB,) ^ . . . . 

Since B is acyclic, Hn{B) is trivial. Therefore, from the exactness of the sequence. 



d : Hn{B,Bc) — i- Hn^i{Bc) is an isomorphism. Let us now use the Theorem 2.4 for 
X = B, A = Be and B = Ba and n = 2- This gives us the isomorphism H2{Ba, Ba n 
Be) = H2{Ba,dBa) ^ Hi{B,Be)- Consequently, we have the sequence of isomorphisms 

H2{Ba,dBa) H2{B,Be) A Hi{Be)- Therefore, the set of cycles . . . , 5(i„}, 

ddi G Ci{Be) for i G {1, . . . , n}, represent a basis of Hi{Be)- 
Now, the non-local Faraday's equations are expressed as 

{t,dd,) = {-zu^,d,),j G {l,...,/3i(/C,)}. (13) 

A novel way to express the jth non local Faraday law in term of the unknowns is to 
pre-multiplying by c-'^j^ 



and using the same passages as when obtaining (11 ) we get 



{iooc^^fiS) n + (c^'^R) T + ^ (c^'^Rc^) = 0. (14) 



^This correspond to algebraically sum local Faraday's equations enforced on dual faces belonging 
to the support of the considered cut. Since the contributions in the interior cancel out, what remains 
is the non local Faraday's law enforced on the boundary of the considered cut, for more details see 
|27j . The cochains on primal complex are denoted by column vectors c*. For the sake of parsimony 
in the notation, the chains di on the dual complex dual to c* are denoted by c'-^. 
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By multiplying the (11) by iuo and considering also the equations (12) and (14) 



the following final symmetric algebraic system having T 
unknowns reads as 



^^\^^,f2andthe {^.j^i^^ 



as 



{iujfi6n + RT + ^^^('^ 



.j=i 



1 



0, 



=0, 



{iujd^^i5) Q, + (c^'^R) T + ^^if (c^'^Rc^) i, 
(T,e) = 0, 



0, jG{l,...,/3i(/CJ}, 
Ve G /Cc \ /Ca- 



(15) 



The source of the problem can be enforced by considering one of the currents ij as 



known, substituting it into (15), and moving its contribution on the right-hand side of 



the system. Alternatively, one can force an e.m.f. by putting its value on the right-hand 
side of the system at the position of the correspondent non-local Faraday's equation, 
see [26l. 



6 A historical survey on the definitions of cuts 

In this Section, three families of definition of cuts presented in the literature in the last 
twenty-five years are reviewed. Due to the use of the Finite Elements with nodal basis 
functions, most definitions concentrate on the so-called thin cuts, which are 2-chains 
on the primal complex. With the modern Finite Elements employing edge elements 
basis functions, cuts defined as in this paper — called thick cuts — are needed in place 
of the thin cuts. Even though it is difficult to generate a thick cut from a thin cut in 
general, see |27], the definitions of thin cuts presented in the following may be easily 
adapted as attempts to define thick cuts also. With nodal basis functions, there is 
the need to impose a potential jump across the thin cuts. This is usually performed 
by "cutting" the cell complex in correspondence of the thin cuts doubling the nodes 
belonging to each cut. This method requires non self-intersecting thin cuts and provides 
some complication when cuts intersect, see for example [53]. On the contrary, the use 
of edge elements, as done in this paper, yields to a straightforward implementation 
even when cuts intersect or, as frequently happens, have self-intersections. 

6.1 Embedded sub- manifolds 

Kotiuga, starting from 1986, published many papers about the definition of cuts, proof 
of their existence and the development of an algorithm to compute them, see 122], |54] . 
[55] . [33] . He defined thin cuts as embedded sub- manifolds being generators of the 
second relative homology group basis H2{fCa,dlCa)- He proposed also an algorithm to 
automatically generate cuts: first a H^{Ba) basis is computed by employing a reduction 
technique based on a tree-cotree decomposition followed by a reordering and a classical 
Smith Normal Form computation [44]. Then, a non-physical Poisson problem is solved 
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for each cut. Finally, cuts are extracted as iso-surfaces of the non-physical problems 
solutions. 

This definition, although a real breakthrough, is too conservative when employing 
the modern edge element basis functions. In fact, in this case there is no need for 
cuts to be embedded surfaces, so the solution of the non-physical problems — which is 
computationally quite costly — can be avoided. 

6.2 Homotopy- based definition 

After Kotiuga's definition, many researchers were persuaded about the existence of 
an easier and more intuitive way to tackle this problem. In our opinion, two reasons 
diverted researchers on heuristic solutions: the first is due to the fact that algebraic 
topology, namely cohomology theory, was — and probably still is — not well known to 
scientists working in computational electromagnetics. The second — perhaps even more 
important — was the lack of efficient algorithms for the computation in a reasonable 
time of cohomology generators. 

In [30], Bott and Tu stated "By some divine justice the homotopy groups or a finite 
polyhedron or a manifold seem as difficult to compute as they are easy to define." In 
fact, a homotopy-based definition of cuts has been introduced becoming soon the most 
popular one. The idea is to introduce a set of 2-cells whose removal transform the 
insulating region into a connected and simply-connected one [56] , [57] , [58] , [59] , |39j^ 
Nonetheless, when dealing with homotopy, one falls easily into intractable problems, 
with a consequent lack of rigorous proofs and details of the algorithms in all the cited 
papers. 

It has been already shown, for example in [60], [61], that there exist cuts that do 
not fulfil the homotopy-based definition. Namely, in case of a knot's complement, the 
cut realization as embedded sub-manifold — which is a Seifert surface — has to leave the 
complement multiply-connected. Even though this counter-example was quite clear, 
dealing with knotted conductors is extremely uncommon in practice. It has been 
concluded that problems in the homotopy-based definition arise only when dealing 
with knot's complement which, as written explicitly by Bossavit [TU p. 238], are really 
marginal in practice. 

Nonetheless, computational electromagnetics community seems not to be aware that 
problems do happen frequently even with the most simple example possible, namely a 
conducting solid torus in which the current fiow. In fact, in this paper we present for 
the first time a concrete counter-example that the homotopy-based definition of cuts 
is not only too restrictive but wrong. What is even more serious is that it makes very 

^This is an attempt for a homotopy-based definition of tliin cuts. Ren in [31 modifies tliis definition 
for tlie tliick cuts. A thick cut is defined as a set of 3-cells whose removal transform the insulating 
region into a simply-connected one. 
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hard even to detect if a potential cut is correct or not, which clearly shows that this 
heuristic definition of cuts should be abandoned. 




Figure 6: The Bing's house counter-example. 

The counter-example is as follows. Consider a solid torus, which represents JCc, and 
its complement /C^ with respect to a ball, which contains the solid torus. By growing an 
acyclic sub-complex and taking the complement, the set of 2-cells in Fig. [6] is obtained. 
Two views are shown in the picture by cutting the set of 2-cells with a vertical plane 
(Fig. [6^) and a horizontal plane (Fig. |6}d). The triangulated torus represents the /Cc 
complex (the torus is not cut with the planes in the picture for the sake of clarity. The 
ball which contains /Cc is not shown either for the same reason.). The set of 2-cells is 
formed by a Bing's house [62\ plus a cylinder. Informally, the torus /Cc is placed in 
the 'upper chamber' of the Bing's house and the hole of the torus is connected to the 
tunnel used to enter the upper chamber by the cylinder. 

It is standard to see that once one removes the set of 2-cells from the complex /C^, 
what remains becomes connected and simply-connected. Hence, this set of 2-cells fulfil 
the homotopy-based definition of cut, but of course this is not a cut since the 1-cycle 
c crosses one time the cut without linking any current. Moreover, it is very hard in 
practice to detect this situation, which makes this definition — and related algorithms — 
not suitable even with the simplest example. 

6.3 Axiomatic definition 

An axiomatic definition of cuts is frequently used in mathematical papers, see for 
example [63] , |6l] , [65] , [66] , [67] , |68] , [69] , [70] , [H^ 

|74j . Cuts are defined as 2-manifolds with boundary {Sj}"^^ which fulfil the following 
set of axioms: 

• The boundary of Sj is located at the boundary of the meshed region in which 
the cuts are searched for. 
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• n = for i ^ j, 

• Q \ ljr=i is pseudo-Lipschitz and simply-connected. 

Once such a set of cuts is provided, the results presented in those papers can be used. 
However, the presented definition of cuts is too restrictive for practical applications. 
One can easily see that even in case of two chained conductors it is not possible to 
find a set of cuts for which Sj fl Sj = holds. The same holds for many practical 
configurations, for example, electric transformers. 

Moreover, this axiomatic definition does not point to an algorithm to find cuts auto- 
matically, which is fundamental for practical problems since it is practically impossible 
to define cuts 'by hand' for serious problems. 

7 Cohomology computation 

A detailed survey on the state-of-the-art algorithms to compute cohomology group gen- 
erators used in electromagnetic modeling can be found in |17]. The solution proposed 
in this paper is to change the available codes for computing homology group generators 
(see |15], |16]) to compute the cohomology group generators. The necessary changes, 
described in detail in |17j, are very easy to implement. 

In order to obtain a computationally efficient code, a so-called shaving procedure for 
cohomology has to be applied. A reduction in (co)homology is a procedure of removing 
from the complex some cells in such a way that the (co) homology groups of the complex 
remains unchanged. Then, the classical Smith Normal Form [44j computation with 
hyper-cubical computational complexity can be performed on the reduced complex. 
A shaving is a reduction of the complex such that the representatives of generators 
in the reduced complex are also representatives of generators in the initial complex. 
As it is explained in [17], the algorithm presented in [75] is a shaving for cohomology 
computations. 

Due to a number of efficient reduction techniques used (namely, [73] followed by 
|76j). in all the tested cases the Smith Normal Form computation has been not used at 
all. The state-of-the-art is the implementation of the acyclic sub-complex shaving with 
look-up tables 07], whose computational complexity is linear and is able to reduce 
almost always the complex down to its cohomology generators. 

An example of cut generated for the complement of a trefoil knot-shaped conductor 
is presented in Fig. [7j 

We would like to point out that it is necessary for the potential design to compute 
the 1st cohomology group generators over integers and this cannot be substituted by 
any field Zp for p prime. In fact, let us consider Z2 as an example. In Fig. [8^, a two 
turn conductor is shown. In Fig. [8]d, the support of a 2-chain dual to a representative 
of a H^{ICa, Z12) generator is presented. When a cycle surrounding the two branches of 
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Figure 7: A trefoil knot conductor together with the dual faces dual to edges belonging 
to the support of the cut. 

the conductor is considered, it does not intersect the support of the chain. It is easy 
to verify that on this cycle the Ampere's law does not hold. Similar examples can be 
constructed for any coefficient field Zp. 




a) b) 

Figure 8: a) A two turn conductor, b) The support of a 2-chain dual to a representative 
of a if^(/Ca5^2) generator. 



8 Conclusions 

In this paper, is has been discussed how the (co)homology theories are fundamental 
for the potential design in computational physics. In particular, a systematic design of 
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potentials employed in the magneto-quasi-static T-Q formulation has been presented. It 
has been demonstrated that the entities called cuts in computational electromagnetics 
are a basis of the first cohomology group over integers of the insulating region. The 
limitations on the definition of cuts presented in the literature are shown by using 
concrete counter-examples, which should persuade the Readers that cohomology is not 
one of the possible options but something which is expressly needed to the potential 
design. 
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